Non-gaussianity from the second-order cosmological perturbation 
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Several conserved and/or gauge invariant quantities described as the second-order curvature per- 
turbation have been given in the literature. We revisit various scenarios for the generation of 
second-order non-gaussianity in the primordial curvature perturbation f, employing for the first 
time a unified notation and focusing on the normalisation /nl of the bispectrum. When £ first 
appears a few Hubble times after horizon exit, |/nl| is much less than 1 and is, therefore, negligible. 
Thereafter £ (and hence /nl) is conserved as long as the pressure is a unique function of energy 
density (adiabatic pressure). Non-adiabatic pressure comes presumably only from the effect of fields, 
other than the one pointing along the inflationary trajectory, which are light during inflation ('light 
non-inflaton fields'). During single-component inflation /nl is constant, but multi-component infla- 
tion might generate |/nl| ~ 1 or bigger. Preheating can affect /nl only in atypical scenarios where 
it involves light non-inflaton fields. The simplest curvaton scenario typically gives /nl -C —1 or 
/nl = +5/4. The inhomogeneous reheating scenario can give a wide range of values for /nl- Unless 
there is a detection, observation can eventually provide a limit |/nl| < 1, at which level it will be 
crucial to calculate the precise observational limit using second order theory. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

Cosmological scales leave the horizon during inflation 
and re-enter it after Big Bang Nucleosynthesis. Through- 
out the super-horizon era it is very useful to define a 
primordial cosmological curvature perturbation, which is 
conserved if and only if pressure throughout the Universe 
is a unique function of energy density (the adiabatic pres- 
sure condition) 0, 0, 0, Gl IS El 01 • Observation directly 
constrains the curvature perturbation at the very end of 
the super-horizon era, a few Hubble times before cosmo- 
logical scales start to enter the horizon, when it appar- 
ently sets the initial condition for the subsequent evolu- 
tion of all cosmological perturbations. The observed cur- 
vature perturbation is almost Gaussian with an almost 
scale-invariant spectrum. 

Cosmological perturbation theory expands the exact 
equations in powers of the perturbations and keeps terms 
only up to the nth order. Since the observed curvature 
perturbation is of order 1CP 5 , one might think that first- 
order perturbation theory will be adequate for all com- 
parisons with observation. That may not be the case 
however, because the PLANCK satellite 0] and its suc- 
cessors may be sensitive to non-gaussianity of the curva- 
ture perturbation at the level of second-order perturba- 
tion theory 0]. 

Several authors have treated the non-gaussianity of 
the primordial curvature perturbation in the context of 
second-order perturbation theory. They have adopted 
different definitions of the curvature perturbation and 



obtained results for a variety of situations. In this pa- 
per we revisit the calculations, using a single definition 
of the curvature perturbation which we denote by £. In 
some cases we disagree with the findings of the original 
authors. 

The outline of this paper is the following: in Section ILTI 
we review two definitions of the curvature perturbation 
found in the literature, which are valid during and af- 
ter inflation, and establish definite relationships between 
them; in section ITTT1 we review a third curvature pertur- 
bation definition, which applies only during inflation, and 
study it in models of inflation of the slow-roll variety; in 
Section llVl we describe the present framework for think- 
ing about the origin and evolution of the curvature per- 
turbation; in Section [V] we see how non-gaussianity is 
defined and constrained by observation; in Section IVTl we 
study the initial non-gaussianity of the curvature pertur- 
bation, a few Hubble times after horizon exit; in Section 
IVIII we study its subsequent evolution according to some 
different models. The conclusions are summarised in Sec- 
tion En] 

We shall denote unperturbed quantities by a subscript 
0, and generally work with conformal time r\ defined by 
the unperturbed line element 



ds 2 = a 2 {rj) (-drj 2 + 5 ij dx l dx : ') 



(1) 



Here a is the scale factor whose present value is taken to 
be 1, and a prime denotes d/di]. Sometimes though we 
revert to physical time t, with a dot meaning d/dt defined 
by d/dt = a~ 1 d/dr). We shall adopt the convention that 
a generic perturbation g is split into a first- and second- 
order part according to the formula 
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II. TWO DEFINITIONS OF THE CURVATURE 
PERTURBATION 

A. Preliminaries 

Cosmological perturbations describe small departures 
of the actual Universe, away from some perfect homo- 
geneous and isotropic universe with the line element 
Eq. ft). For a generic perturbation it is convenient to 
make the Fourier expansion 

5(x, V) = j^y^ I d 3 k g(k, V y k * , (3) 

where the spacetime coordinates are those of the unper- 
turbed Universe. The inverse of the comoving wavenum- 
ber, fc _1 , is often referred to as the scale. 

A given scale is said to be outside the horizon during 
the era aH 3> k 7 where H = a/a is the Hubble parame- 
ter. Except where otherwise stated, our discussion applies 
only to this super-horizon regime. 

When evaluating an observable quantity only a lim- 
ited range of scales will be involved. The largest scale, 
relevant for the low multipoles of the Cosmic Microwave 
Background anisotropy, is fc _1 ~ Hq 1 where Hq is 
the present Hubble parameter. The smallest scale usu- 
ally considered is the one enclosing matter with mass 
~ 10 6 M©, which corresponds to fc _1 ~ 10~ 2 Mpc ~ 
IO^Hq 1 . The cosmological range of scales therefore ex- 
tends over only six decades or so. 

To define cosmological perturbations in general, one 
has to introduce in the perturbed Universe a coordinate 
system (t,x l ), which defines a slicing of spacetime (fixed 
t) and a threading (fixed x l ). To define the curvature 
perturbation it is enough to define the slicing |?J . 

B. Two definitions of the curvature perturbation 

In this paper we take as our definition of C the follo wing 
expression for the spatial metric |1 U. IToL ill 111 HI Il4f 

which applies non-perturbatively: 

9ij = a 2 (?7)7ye 2C . (4) 

Here jij has unit determinant, and the time-slicing is one 
of uniform energy density 1 . 

It has been shown under weak assumptions Q that 
this defines C uniquely, and that ( is conserved as long 
as the pressure is a unique function of energy density. 
Also, it has been shown that the uniform density slicing 



1 It is proved in Ref. Q tha t this definition of £ coincides with 
that of Lyth and Wands ll , provided that their slices of uniform 
coordinate expansion are taken to correspond to those on which 
the line element has the form Eq. J1J without the factor e 2 ^ (this 
makes the slices practically flat if 7ij ~ 



practically coincides with the comoving slicing (orthogo- 
nal to the flow of energy) , and with the uniform Hubble 
slicing (corresponding to uniform proper expansion, that 
expansion being practically independent of the threading 
which defines it) The coincidence of these slicings is 
important since all three have been invoked by different 
authors. 

Since the matrix 7 has unit determinant it can be writ- 
ten t_ = Ie h , where / is the unit matrix and h is trace- 
less |3- Assuming that the initial condition is set by 
inflation, h corresponds to a tensor perturbation (gravi- 
tational wave amplitude) which will be negligible unless 
the scale of inflation is very high. As we shall see later 
(see footnote j the results we are going to present are 
valid even if h is not negligible, but to simplify the pre- 
sentation we drop h from the equations. Accordingly, the 
space part of the metric in the super-horizon regime is 
supposed to be well approximated by 

9ij = a 2 (?y)<5y e 2C . (5) 

At first order, Eq. (0 corresponds to 

ft, =a 2 (7 7 )%(l + 2C). ( 6 ) 

Up to a sign, this is the definition of the first-order cur- 
vature perturbation adopted by all authors. There is no 
universally agreed convention for the sign of Ours 
coincides with the convention of most of the papers to 
which we refer, and we have checked carefully that the 
signs in our own set of equations are correct. 
At second order we have 

m =a 2 (77)Ml + 2C + 2C 2 ). (7) 

This is our definition of £ at second order. 

Malik and Wands [0| instead defined ( by Eq. © even 
at second order. Denoting their definition by a subscript 
MW, 

C MW = C + C 2 , (8) 

or equivalently 

C 2 MW = C2+2(Ci) 2 , (9) 

where £1 is the first-order quantity whose definition 
Eq. 10 is agreed by all authors. 

To make contact with calculations of the curvature per- 
turbation during inflation, we need some gauge-invariant 
expressions for the curvature perturbation. 'Gauge- 
invariant' means that the definition is valid for any choice 
of the coordinate system which defines the slicing and 
threading 2 . 



2 In the unperturbed limit the slicing has to be the one on which all 
quantities are uniform and the the threading has to be orthogonal 
to it. 
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We shall write gauge-invariant expressions in terms of 
£ and £ MW . First we consider a quantity ip MW , denned 
even at second order by 



g ij = c?(r } )5 ij {l-2^) 



(10) 



This definition, which is written in analogy to Eq. Jfjjl, 
applies to a generic slicing. Analogously to Eq. JSJ we 
can consider a quantity tp, valid also in a generic slicing, 
defined by 



III. SLOW-ROLL INFLATION AND A THIRD 
DEFINITION 

Now we specialize to the era of slow-roll inflation. 
We consider single-component inflation, during which 
the curvature perturbation £ is conserved, and multi- 
component inflation during which it varies. After defin- 
ing both paradigms, we give a third definition of the cur- 
vature perturbation which applies only during inflation. 



9ij = a 2 (r))5ije' 



-2ip 



(11) 



On uniform-density slices, ip t = = V^ W = 

-( 2 MW , and ip2 — —(,2- We shall also need the energy 
density perturbation 5p, defined on the generic slicing, 
as well as the unperturbed energy density po. 

At first order, the gauge-invariant expression for £ has 
the well-known form 



Ci 



-ipi - H— 
Po 



(12) 



where Ji = a'/a, and the unperturbed energy density 
satisfies /3q = — 37i(po + -Po) with Po being the unper- 
turbed pressure. This expression obviously is correct for 
the uniform density slicing, and it is correct for all slic- 
ings because the changes in the first and second terms 
induced by a change in the slicing cancel 0, 0, 0, E3j E3 • 
At second order, Malik and Wands show that Hfl 
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Po 



2H 
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(13) 



which is, again and for the same reason as before, obvi- 
ously correct for all the slices. Accordingly, from Eq. 
we can write a gauge invariant definition for our second- 
order £: 3 



(2 
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where the relation 



V 2 MW -^2-2(^) 2 , 
coming from Eqs. I|10[) and i|llfl. has been used. 



(14) 



(15) 



3 This relation has recently been confirmed in Ref. 
nonlinear coordinate-free approach. 



A. Single-component inflation 



In a single-component inflation model |2(J, |21j the in- 
flaton trajectory is by definition essentially unique. The 
inflaton field tp parameterises the distance along the in- 
flaton trajectory. In terms of the field variation, slow-roll 
inflation is characterised by the slow-roll conditions e <C 1 
and \r] V \ < 1 US El, where 



Vv 



= -H/H 2 , 

_ <£_ 

Hip ' 



(16) 
(17) 



The inflaton field can be taken to be canonically nor- 
malised, in which case these definitions are equivalent to 
conditions on the potential V 



M 2 



8VY 



2V 2 \dp J 
M 2 d 2 V 

V dip 2 



(18) 
(19) 



which, together with the slow-roll approximation, lead to 
the slow-roll behaviour 



3Hip 



dV 

dip 



(20) 



Here Mp is the reduced Planck mass (Mp = (SttGn)^ 1 ). 

Even without the slow-roll approximation, slices of uni- 
form tp correspond to comoving slices because a spatial 
gradient of tp would give non-vanishing momentum den- 
sity. Since comoving slices coincide with slices of uniform 
energy density, the slices of uniform tp coincide also with 
the latter. Also, since tp is a Lorentz scalar, its gauge 
transformation is the same as that of p. It follows |23 
that we can replace p by tp in the above expressions: 



Ci 



S>2 



-Ipl - U—f- 
Vo 



(21) 
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Po 
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(22) 
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~4>2 



H 



go 



(23) 



B. Multi-component inflation 

Now consider the case of multi-component inflation, 
where there is a family of incquivalcnt inflationary trajec- 
tories lying in an TV-dimensional manifold of field space. 
If the relevant part of the manifold is not too big it will 
be a good approximation to take the fields to be canoni- 
cally normalised. Then the inequivalent trajectories will 
be curved in field space 4 . To define the trajectories one 
can choose a fixed basis in field space corresponding to 
fields 0i, • ■ • , 4>n- 

Assuming canonical normalisation, multi-component 
slow-roll inflation is characterised by the conditions 



Ml 
2V 2 
Ml 
V 



dV 



d 2 V 



•C 



(24) 

(25) 
(26) 




const 



= const 



V = const 



(a) 




const 



V = const 



V = const 



FIG. 1: Two different procedures for defining the fields in 
two-component inflation. The fields are denoted by tp and a. 
(a) The field tp parameterises the distance along the inflaton 
trajectories, with uniform tp corresponding to the equipoten- 
tial lines. The field a parameterises the distance along the 
equipotentials. (b) The fields tp and a are the components in 
a fixed orthonormal basis, aligned with the inflationary tra- 
jectory at a certain point in field space. The value of tp is now 
the displacement along the tangent vector and the value of 
a is the displacement along the orthogonal vector. Working 
to second order in these displacements, the equipotentials no 
longer coincide with the lines of uniform tp. 



The procedure of choosing a fixed basis is quite conve- 
nient for calculations, but a different procedure leads to a 
perhaps simpler theoretical description. This is to take (p 
to parameterise the distance along the inflaton trajecto- 
ries, just as in single-component inflation, but now with 
the proviso that uniform tp corresponds to uniform field 
potential (since we work in the slow-roll approximation, 
this means that the slices in field space of uniform (p are 
orthogonal to the trajectories). Then, in the slow-roll 
approximation, slices of spacetime with uniform tp will 
again coincide with slices of uniform density (see Fig. 
Oil). Since tp is a scalar, Eqs. I|21|) and (|22|l will then be 
valid. This is the simplest form of the gauge-invariant 
expression, though for a practical calculation it may be 
better to write it in terms of a fixed basis. 

There is a subtlety here. For the first-order case we 
could define tp in a different way; around a given point 
on the unperturbed trajectory we could choose a fixed 
field basis, with one of the basis vectors pointing along 
the trajectory, and define tp as the corresponding field 
component. Then we could choose tp to be canonically 
normalised in the vicinity of the chosen point in field 
space. That would not work at second order though, 
because at that order it makes a difference whether tp is 
the appropriate parameterisation of the distance along 



4 More generally they will be non-geodesics, the geodesies being 
the trajectories which the background fields could follow if there 
was no potential term in the scalar Lagrangian l23l . 



the trajectories (our adopted definition) or the distance 
along a tangent vector to the trajectory (the alternative 
definition) (see Fig. ^p) . Only our adopted one will make 
Eqs. (|22Jl and (J23J) valid. 



C. A third definition of the curvature perturbation 

The third definition in the literature applies only dur- 
ing inflation. It was given originally by Acquaviva et. al. 
|24j for the single-component case, and the generalization 
to the multi-component case was noted by Rigopoulos 
|23|. We shall denote this definition by £ A . 

The definition of Acquaviva et. al. and Rigopoulos is 



C 2 A = 



(^ + 2^1 + WiM) 

W + 2H 2 - HtpZ/tp'v 



(27) 



This is gauge-invariant by construction, with tp defined 
as in Figure 1(a). 

It was pointed out by Vernizzi |22| (actually in the con- 
text of single-component inflation) that comparing this 
definition with Eq. 1)13(1 gives simply 



C 2 A 



S>2 



W 2 (Ci) 2 



W + 2H 2 - Htp"/tp' 



(28) 
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In the limit of slow-roll the denominator of the last term 
becomes just 2TL 2 , and then 

C 2 A = C2. (29) 

In other words, this third definition coincides with our 
adopted one in the slow-roll limit. 

Making use of the slow-roll parameters defined in Eqs. 
(|16fl and (|17fl . the expression in Eq. I|28|l gives to first- 
order in the slow-roll approximation 

C 2 A = C2-(2e-w)(Ci) 2 . (30) 

IV. THE EVOLUTION OF THE CURVATURE 
PERTURBATION 

The simplest possibility for the origin of the observed 
curvature perturbation is that it comes from the vacuum 
fluctuation of the inflaton field in a single-component 
model. More recently other possibilities were recognised 
and we summarise the situation now. Although the 
discussion is usually applied to the magnitude of the 
curvature perturbation, it applies equally to the non- 
gaussianity. 

A. Heavy, light and ultra-light fields 

On each scale the initial epoch, as far as classical per- 
turbations are concerned, should be taken to be a few 
Hubble times after horizon exit during inflation. The 
reason is that all such perturbations are supposed to orig- 
inate from the vacuum fluctuation of one or more light 
scalar fields, the fluctuation on each scale being promoted 
to a classical perturbation around the time of horizon 
exit. 

Considering a fixed basis with canonical normalisation, 
a light field is roughly speaking one satisfying the flat- 
ness condition in Eq. 1251) . The terminology is suggested 
by the important special case that the effective poten- 
tial during inflation is quadratic. Then, a light field is 
roughly speaking that whose effective mass during infla- 
tion is less than the value H„ of the Hubble parameter. 
More precisely, the condition that the vacuum fluctuation 
be promoted to a classical perturbation is [2j| 

3 

m < -H* . (31) 

From now on we focus on the quadratic potential, and 
take this as the the definition of a light field. Conversely a 
heavy field may be defined as one for which the condition 
in Eq. (|3"TJl is violated. 

During inflation light fields slowly roll according to 
Eq. H26fl (with the vacuum fluctuation superimposed) 
while the heavy fields presumably are pinned down at 
an instantaneous minimum of the effective potential. As 
we have seen, multi-component inflation takes place in a 



subspace of field space. The fields in this subspace are 
light, but their effective masses are sufficient to apprecia- 
bly curve the inflationary trajectories. In the case of both 
multi-component and single-component inflation, there 
could also be 'ultra-light' fields, which do not apprecia- 
bly curve the inflationary trajectory and which therefore 
have practically no effect on the dynamics of inflation. 

B. The evolution of the curvature perturbation 

To describe the behaviour of perturbations during the 
super-horizon era, without making too many detailed as- 
sumptions, it is useful to invoke the separate universe hy- 
pothesis 0, IE 13 151 after smoothing on a given comov- 
ing scale much bigger than the horizon 5 . According to 
this hypothesis the local evolution at each position is that 
of some unperturbed universe (separate universe). Of 
course the separate universe hypothesis can and should 
be checked where there is a sufficiently detailed model. 
However, it should be correct on cosmological scales for 
a very simple reason. The unperturbed Universe may 
be defined as the one around us, smoothed on a scale 
a bit bigger than the present Hubble distance. In other 
words, the separate universe hypothesis is certainly valid 
when applied to that scale. But the whole range of cos- 
mological scales spans only a few decades. This means 
that cosmological scales are likely to be huge compared 
with any scale that is relevant in the early Universe, and 
accordingly that the separate universe hypothesis should 
be valid when applied to cosmological scales even though 
it might fail on much smaller scales(this expectation was 
verified in a preheating example |27| to which we return 
later). 

We are concerned with the curvature perturbation, 
which during the super-horizon era is conserved as long 
as the pressure is a unique function of the energy density 
(the adiabatic pressure condition). The adiabatic pres- 
sure condition will be satisfied if and only if the separate 
universes are identical (at least as far as the relation be- 
tween pressure and energy density is concerned) 6 . The 
condition to have identical universes after a given epoch 
is that the specification of a single quantity at that epoch 
is sufficient to determine the entire subsequent evolution. 

In the case of single-component inflation, the initial 
condition may be supplied by the local value of the in- 
flaton field, at the very beginning of the super-horizon 



5 When considering linear equations, smoothing is equivalent to 
dropping short wavelengths fourier components. In the nonlin- 
ear situation the smoothing procedure could be in principle am- 
biguous. In a given situation one should state explicitly which 
quantities are being smoothed. 

6 Of course the identity will only hold after making an appropriate 
synchronization of the clocks at different positions. Having made 
that synchronization, horizon entry will occur at different times 
in different positions, which can be regarded as the origin of the 
curvature perturbation. 
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era when it first becomes classical. Given the separate 
universe hypothesis, that is the only possibility if the in- 
flaton is the only light field ever to play a significant dy- 
namical role. This means that the curvature perturbation 
generated at horizon exit during single- component infla- 
tion will be equal to the one observed at the approach of 
horizon entry, provided that the inflaton is the only light 
field ever to play a dynamical role. 

If inflation is multi-component, more than one field is 
by definition relevant during inflation. Then the curva- 
ture perturbation cannot be conserved during inflation. 
The variation of the curvature perturbation during multi- 
component inflation is caused by the vacuum fluctuation 
orthogonal to the unperturbed inflationary trajectory, 
which around the time of horizon exit kicks the trajectory 
onto a nearby one so that the local trajectory becomes 
position-dependent. After inflation is over, the curvature 
perturbation will be conserved if the local trajectories 
lead to practically identical universes. In other words it 
will be conserved if the light (and ultra-light) fields, or- 
thogonal to the trajectory at the end of inflation, do not 
affect the subsequent evolution of the Universe. 

The curvature perturbation after inflation will vary if 
some light or ultra-light field, orthogonal to the trajec- 
tory at the end of inflation, affects the subsequent evolu- 
tion of the Universe (to be precise, affects the pressure). 
As we shall describe in Section lVlII three types of scenar- 
ios have been proposed for this post-inflationary variation 
of the curvature perturbation. 



V. NON-GAUSSIANITY 

A. Defining the non-gaussianity 

A gaussian perturbation is one whose Fourier compo- 
nents are uncorrelated. All of its statistical properties 
are defined by its spectrum, and the spectrum Vn of a 
generic perturbation is conveniently defined pol . |21| by 7 

( 5 (ki)g(k 2 )) = ^ 3 (kx + k 2 )^(fc) , (32) 
the normalisation being chosen so that 



<5 2 (x)> 



dk 



(33) 



On cosmological scales a few Hubble times before horizon 
entry, observation shows that the curvature perturbation 
is almost Gaussian with V^ 2 ~ 10~ 5 . 



7 Technically the expectation values in this and the following ex- 
pressions refer to an ensemble of universes but, because the 
stochastic properties of the perturbations are supposed to be 
invariant under translations, the expectation values can also be 
regarded as averages over the location of the observer who defines 
the origin of coordinates. 



The simplest kind of non-gaussianity that the curva- 
ture perturbation could possess is of the form 



C(x) 



Cg(x) - -/nl 



(34) 



where ( g is Gaussian with ((g) = 0, and the non-linearity 
parameter /nl is independent of position. We will call 
this correlated \ 2 non-gaussianity. Note that this defi- 
nition assumes that (() =0, which means that the zero 
Fourier mode (spatial average) is dropped. 

Following Maldacena ^(|> we have inserted the pref- 
actor —(3/5) so that in first-order perturbation theory 
our definition agrees with that of Komatsu and Spergel 
0, which is generally the definition people use when 
comparing theory with observation. Working in first- 
order perturbation theory, these authors write ^(x) = 

<I> g (x) + /nl ( $g(x) - <&j*J 5 an d their $ is equal to -3/5 

times our (. 

One of the most powerful observational signatures of 
non-gaussianity is a nonzero value for the three-point cor- 
relator, specified by the bispectrum B defined by pj.l2ftj| 

(C(k 1 )C(k 2 )C(k 3 )) = (2TT)- 3 / 2 B(k 1 ,k 2 ,k 3 )6 3 (k 1 +k 2 +k 3 ). 

For correlated \ 2 non-gaussianity (with the gaussian 
term dominating) 



B(k!,k 2 ,k 3 ) = --/nl 
5 



P c (ki)P c (k 2 ) + cyclic , (36) 



where -P<j(fc) = 2ir 2 'P ( ;(k)/k 3 . For any kind of non- 
Gaussianity one may use the above expression to define 
a function /nlCsi, k 2 , k 3 ). 

Given a calculation of /nl using first-order perturba- 
tion theory, one expects in general that going to second 
order will change /nl by an amount of order I. On this 
basis, one expects that a first-order calculation is good 
enough if it yields |/nl| 3> 1, but that otherwise a second- 
order calculation will be necessary. 

The definition Eq. (|36fl of /nl is made using our 
adopted definition of (. If Q in the definition is replaced 
by £ MW (with the zero Fourier mode dropped) then /nl 
should be replaced by 



f MW _ t 
/NL = /NL - - • 



(37) 



To obtain this expression we used Eq. ijHJ and dropped 
terms higher than second order. 

All of this assumes that the non-gaussian component 
of C is fully correlated with the gaussian component. An 
alternative possibility |29j that will be important for us 
is if C has the form 



C(x)=c s (x)-|/ NL (C(x)-c, 



>2 ' 



(38) 



where ( 9 and ( CT are uncorrelated Gaussian perturbations, 
normalised to have equal spectra, and the parameter /nl 
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is independent of position. We will call this uncorrelated 
\ 2 non-gaussianity. It can be shown |29j that in this 
case, /nl as defined by Eq. ijjffijl is given by 



parameter e: 



/NL 



XL 



1300 



(39) 



B. 



Observational constraints on the 
non-gaussianity 



Taking /nl to denote the non-linearity parameter at 
the primordial era, let us consider the observational con- 
straints. Detailed calculations have so far been made only 
with /nl independent of the wavenumbers, and only by 
using first-order perturbation theory for the evolution of 
the cosmological perturbations after horizon entry. It is 
found 30] that present observation requires |/nl| < 10 2 
making the non-gaussian fraction at most of order 10~ 3 . 
The use of first-order perturbation theory in this con- 
text is amply justified. Looking to the future though, 
it is found that the PLANCK satellite will either detect 
non-gaussianity or reduce the bound to |/nl| < 5 HIHI, 
and that foreseeable future observations can reach a level 

i/nl| ~3 UM 

Although the use of first-order perturbation theory is 
not really justified for the latter estimates, we can safely 
conclude that it will be difficult for observation ever to de- 
tect a value |/nl| <SC 1. That is a pity because, as we shall 
see, such a value is predicted by some theoretical scenar- 
ios. On the other hand, other scenarios predict |/nl| 
roughly of order 1. It will therefore be of great interest 
to have detailed second-order calculations, to establish 
precisely the level of sensitivity that can be achieved by 
future observations. A step in this direction has been 
taken in Refs. [3ll l32| . where the large-scale cosmic mi- 
crowave background anisotropy is calculated to second 
order in terms of only the curvature perturbation (gen- 
eralizing the first-order Sachs- Wolfe effect H^|). 



VI. THE INITIAL NON-GAUSSIANITY 

A. Single-component inflation 

At first order, the curvature perturbation during 
single-component inflation is Gaussian. Its time- 
independent spectrum is given by [2(], 0] 



P c (k) = 



H\ 2 ( H x :1 



2tt 



(40) 



J k=aH 

and its spectral index n = d In Vq /d In k is given by 

n - 1 = 2r? y - 6e . (41) 

The spectrum r of the tensor perturbation, defined as 
a fraction of Vq 1 is also given in terms of the slow-roll 



r = 16e. 



(42) 



// the curvature perturbation does not evolve after 
single- component inflation is over observation constrains 
n and r, and hence the slow-roll parameters r\y and e. 
A current bound is -0.048 < n - 1 < 0.016 and 
r < 0.46. The second bound gives e < 0.029, but barring 
an accurate cancellation the first bound gives e < 0.003. 
In most inflation models e is completely negligible and 
then the first bound gives —0.024 < r\ v < 0.008 (irre- 
spective of slow-roll inflation models, the upper bound in 
this expression holds generally, and the lower bound is 
badly violated only if there is an accurate cancellation) . 
The bottom line of all this is that e and \r]v\ are both 
constrained to be < 10~ 2 . 

Going to second order, Maldacena 01 nas calculated 
the bispectrum during single-component inflation (see 
also Refs. Elll III III HI IH HI). His result may be 
written in the form 



/nl = (2w 



6e-2ef(k 1 ,k 2 ,k 3 )) , (43) 



with < / < 5/6. By virtue of the slow-roll conditions, 
/nl| "C 1 8 . In other words, the curvature perturbation 
£, as we have defined it, is almost Gaussian during single- 
component inflation. 

From Eq. (|30|l £ A is also practically gaussian, but this 
quantity is defined only during inflation and therefore 
could not be considered as a replacement for £. More im- 
portantly, £ MW has significant non-gaussianity because, 
from Eq. it corresponds to « -5/3. 

One may ask why it is our £ and not £ MW which is 
gaussian in the slow-roll limit 9 . One feature that distin- 
guishes our is that any part of it can be absorbed into 
the scale factor without altering the rest; indeed 



9ij = S lj a 2 {n)e 



2C1+C2 



S tj a (rj)e 



C2 



(44) 



with a = ae 1 * 1 (if we tried to do that with £mWj the 
part of C not absorbed would have to be re-scaled) . This 
means that an extremely long-wavelength and possible 
large part of £ has no local significance. It also means, 
in the context of perturbation theory, that the first-order 
part of £ can be absorbed into the scale factor when dis- 
cussing the second-order part. However, the gaussianity 
of ( does not seem to be related directly to this feature. 



8 Near a maximum of the potential 'fast-roll' inflation |3flt l4flj| can 
take place with \r]y \ somewhat bigger than 1. Maldacena's cal- 
culation does not apply to that case but, presumably, it gives 
initial non-gaussianity |/nl| ~ 1- However, the precise initial 
value of /nl in this case is not important because the corre- 
sponding initial spectral index is far from 1, which means that 
the initial curvature perturbation must be negligible. 

9 We thank Paolo Creminelli for enlightening correspondence on 
this question. 
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Rather, it has to do with the gauge transformation, re- 
lating quantities ipA and ips defined on different slicings. 
With our definition Q, the gauge transformation is 

^(t,x)-V B (t,x) = -AN AB (t,x), (45) 

where AN A b is the number of e-folds of expansion going 
from a slice B to a slice A, both of them corresponding to 
time t 10 . In writing this expression we used physical time 
t instead of conformal time, the two related by dt = adr\. 
Along a comoving worldline, the number of e-folds of 
expansion is defined as N = J Hdr where H is the local 
Hubble parameter and dr is the proper time interval Q • 
To understand the relevance of this result, take ipB = 
and ipA = — C The pressure is adiabatic during single- 
component inflation, which means that dt can be iden- 
tified with the proper time interval dr. and the proper 
expansion rate on slicing A is uniform Q. As a result, to 
second order, 

C = H(t)At(t,x) + ^H(t)(At(t,x)) 2 

~ HAt(t,x) + ^(HAt(t,x)) 2 

~ HAt(t,x). (46) 

In the last line we made the slow-roll approximation, and 
from the second line we can see that the error in /nl 
caused by this approximation is precisely e. 

We also need the gauge transformation for the inflaton 
field ip in terms of At. Since the slices correspond to the 
same coordinate time, the unperturbed inflaton field can 
be taken to be the same on each of them which means 
that the gauge transformation for Sip is 

5ip A (t,x) -6<p B (t,x) =Aip AB (t,x), (47) 

where Aip A B is the change in ip going from slice B to slice 
A. But slice A corresponds to uniform ip., which means 
that on slice B to second order 

ff(t)^ = -H(t)At{^)-\H{t)^{At{t^)f 

<Po 2 ip 

~ -HAt{t^)- l --^-{HAt{t,x)f 
2 Hip 

~ -HAt(t,x), (48) 

where in the last line we used the slow-roll approxima- 
tion. We can see that the fractional error caused by this 
approximation is ipo/Hipo — e — r\y . 

Combining Eqs. (|46|l and i|48|) we have in the slow-roll 
approximation 

t«-ir(t)^, (49) 



This expression is valid even when the tensor perturbation is 
included 0. As a result, the gauge-invariant expressions men- 
tioned earlier are still valid in that case, as are the results based 
on them including the present discussion. 



with fractional error of order max{ ryy,e} (this can also 
be seen directly from Eqs. (|21|1 and H23[l evaluated with 
■0 = 0, but we give the above argument because it ex- 
plains why the result is valid for £ as opposed to £ MW ). 

The final and crucial step is to observe that in the slow- 
roll approximation ip B is gaussian, with again a fractional 
error of order max{?7y,e}. This was demonstrated by 
Maldacena [10( but the basic reason is very simple. The 
non-gaussianity of ip comes either from third and higher 
derivatives of V (through the field equation in unper- 
turbed spacetime) or else through the back-reaction (the 
perturbation of spacetime); but the first effect is small 
[2fi l21j by virtue of the flatness requirements on the po- 
tential, and the second effect is small because (po/H 2 is 
small 21]. This explains why £ with our adopted defi- 
nition is practically Gaussian by virtue of the slow-roll 
approximation. 



B. Multi-component inflation 

The flatness and slow-roll conditions Eqs. Q24[). I|25|l . 
and J2U ensure that the curvature of the inflationary 
trajectories is small during the few Hubble times around 
horizon exit, during which the quantum fluctuation is 
promoted to a classical perturbation. As a result, the 
initial curvature perturbation in first-order perturbation 
theory, is still given by Eq. H4U|) in terms of the field ip 
that we defined earlier. 

What about the initial non-gaussianity generated at 
second order? In the approximation that the curvature 
of the trajectories around horizon exit is completely neg- 
ligible, the orthogonal fields are strictly massless. Such 
fields would not affect Maldacena's second-order calcu- 
lation, which would therefore still give the initial non- 
gaussianity. It is not quite clear whether the curvature 
can really be neglected, as it can be in the first-order case, 
and it may therefore be that the initial non-gaussianity 
in multi-field inflation is different from Maldacena's re- 
sult. Even if that is the case though, we can safely say 
that the initial non-gaussianity corresponds to |/nl| <C 1, 
since the curvature of the trajectories is certainly small. 

VII. THE EVOLUTION AFTER HORIZON EXIT 
A. Single-component inflation and <^ 

During single-component inflation the curvature per- 
turbation C, as we have defined it, does not evolve. From 
its definition Eq. (JSJ, the same is true of C2 ■ 

In contrast C^S given by Eq. (|30(l . will have the slow 
variation 

4 A «-(26-w)(d) 2 . (50) 

This variation has no physical significance, being an ar- 
tifact of the definition. 
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Using a particular gauge, Acquaviva et. al. |24( have 
calculated (2 m terms of first-order quantities ipi, 8<p\, 
and their derivatives, and they have displayed the result 
as an indefinite integral 



cHt) 



A(t)dt + B{t) . 



(51) 



Inserting an initial condition, valid a few Hubble times 
after horizon exit, this becomes 



C 2 A W = CHU) + / A{t)dt+B\\ 



(52) 



In view of our discussion, it is clear that these equations 
will, if correctly evaluated, just reproduce the time de- 
pendence of Eq. (|5U|l . 

The authors of Ref. [24| also present an equation for 
£2 , again involving only first-order quantities, which is 
valid also before horizon entry. Contrary to the claim 
of the authors, this classical equation cannot by itself be 
used to calculate the initial value (more precisely, the 
stochastic properties of the initial value) of (2- In par- 
ticular, it cannot by itself reproduce Maldacena's calcu- 
lation of the bispectrum. 

It is true of course that in the Heisenberg picture the 
quantum operators satisfy the classical field equations. 
In first-order perturbation theory, where the equations 
are linear, this allows one to calculate the curvature per- 
turbation without g oing to the trouble of calculating the 
second-order action |2l| (at the nth order of perturbation 
theory the action has to be evaluated to order n + 1 if it 
is to be used) . At second order in perturbation theory it 
remains to be seen whether the Heisenberg picture can 
provide a useful alternative to Maldacena's calculation, 
who adopted the interaction picture and calculated the 
action to third order. 



B. Multi-component inflation 

During multi-component inflation the curvature per- 
turbation by definition varies significantly along a generic 
trajectory, which means that non-gaussianity is gener- 
ated at some level. So far only a limited range of models 
has been investigated [H El El El El El . To keep 
the spectral tilt within observational bounds, the unper- 
turbed trajectory in these models has to be specially cho- 
sen, but the choice might be justified by a suitable initial 
condition. 

We shall consider here a calculation by Enqvist and 
Valhkoncn in Ref. |45| . Following the same line as Ac- 
quaviva et. al. [24(, they study a two-component in- 
flation model, in which the only important parts of the 
potential are 



V{ip,<r) =V + 



2 2 
m„cr 



1 



2 2 



(53) 



The masses are both supposed to be less than (3/2)ii», 
so that this is a two-component inflation model, and the 
above form of the potential is supposed to hold for some 
number AN of e-folds after cosmological scales leave the 
horizon. They take the unperturbed inflation trajec- 
tory to have Co = 0, and the idea is to calculate the 
amount of non-gaussianity generated after AiV e-folds. 
Irrespective of any later evolution, this calculated non- 
gaussianity will represent the minimal observed one (un- 
less non-gaussianity generated later happens to cancel 
it). 

It is supposed that the condition er = 0, as well as 
the ending of inflation, will come from a tree-level hybrid 
potential, 



V(<p,a) =V - \ml<r 2 + ^Xa 4 



1 2 2 



1 2 2 2 
2 9(Jlp 



,( 54 ) 

Like the original authors though, we shall not investigate 
the extent to which Eq. I|54|l can reproduce Eq. I|53[) for at 
least some number of e-folds. We just focus on Eq. Q53p. 
with the assumption 00 = for the unperturbed trajec- 
tory. 

Because <tq = 0, the unperturbed trajectory is straight, 
and at first order the curvature perturbation £ is con- 
served. This is not the case though at second order. 
Adopting the definition £ A , the authors of El] give an 
expression for (2 similar to that in Eq. (|52ll describing 
the evolution of the second-order curvature perturbation 
on superhorizon scales 11 . This equation, in the gener- 
alized longitudinal gauge, reads (from Eq. (67) in Ref. 
I!): 



C 2 A W-C 2 A (ti) 



1 



eHM 2 



+4A~ 2 d l {6& 1 d l Sa 1 y +ml{Sa 1 ) 2 
+(e - r ]v )6HA- i d l (d k d k 8a 1 d l S<j 1 y 
+ (e - rj v )HA- 4 d i d i (d k 8aid k 8aiy 
-3A-%{d k d k 5aid i Sa 1 y 



~A- i d i d i {d k 5a 1 d k 8a 1 y 



dt 



- A~ 2 d i (S<Tid i 5(Ji) + 3A- 4 d i (d k d k 8a 1 d i 5<Tiy 
-A- 4 d l d l {d k 8a 1 d k 8(T 1 y 



-3eHA-%{d k d k 8a 1 d t Sa 1 



eH 



%d i (d k Sa 1 a k Strij\\ t t } , 



(55) 



11 The fields ip and u in Eq. 1531 are supposed to be canonically 
normalised, which means that ip is not the field appearing in 
th e R igopoulos definition Eq. 1271 of f A . Instead the authors 
of 45J give an equivalent definition in terms of the canonically 
normalised fields. 
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where A~ 2 is the inverse of the Laplacian operator. 

Assuming that this expression is correct, we consider 
the non-gaussianity it may generate. Following the orig- 
inal authors, we note first that 



ICi 



constant, 



(56) 



which is a good approximation since the first-order per- 
turbation equation of the effectively massless field a is 
the same as the first-order perturbation equation of the 
inflaton field cp on superhorizon scales. Moreover, the 
time derivative of the first-order perturbation in a can 
be estimated as 



\6&x\ 



H 



\Sa± 



(57) 



assuming slow-roll conditions. If we also assume that H 
and m 2 are almost constants in time, we end up with 



C 2 A W-C 2 A (*0 = - 



1 



eHM 2 j i 



6HAT 3 dt(5&i8 i 6ai) 
2{6<7 1 ) 2 + ml{Sa 1 ) 2 ]dt. (58) 



For m 2 / H 2 <C 1 the first and third term in the inte- 
grand dominate, whereas all the three terms become of 
the same order of magnitude for m 2 /H 2 ~ 1 which is the 
limit of applicability of the calculation. In any case, the 
typical magnitude of the right-hand side is of order 



-2(fci) 2 + mJ(foi) : 



dt 



eHM 2 Jt 



-2(Scnf 



2A- 2 d i (6a 1 d l 6a 1 ) 
(60) 



where the equation of motion <5<7i + 3HSai + m 2 8ai = 
has been used. Enqvist and Valhkonen seem to have 
neglected the factor 8&\ in the above expression, keeping 
only the term —2(8&i) 2 in the integrand, which gives the 
partial estimate 



C A (i)-C A (*i) 



<< 



H 2 



ICi 



(61) 



This last step is wrong because 8d\ can only be neglected 
compared with 3HS&\ + mr a 8a\, which is not present in 
the integrand in Eq. H60|) . Moreover, what we need is 
the integral of 8a\, i.e. 8&i, which is in any case non- 
negligible. That is the origin of our new term 



AiV^ICrl 



(62) 



in Eq. (|59|l . As a conclusion the level of non-gaussianity 
in this hybrid-type model seems to be much smaller than 
previously thought as the condition m a < H is presum- 
ably not satisfied for a sufficient number of e-folds. 



C. Preheating 



A^lCrl 2 



(59) 



with AN the number of e-folds specified by the integral 
limits. This looks big, but we have to remember that the 
right hand side is uncorrelated with the inflaton pertur- 
bation 8(f) which generates Therefore, Eq. I)38fl as op- 
posed to Eq. H34f> applies, and we would need AN ~ 1300 
to get even /nl ~ 1, which is impossible. 

These conclusions differ sharply from those of En- 
qvist and Vaihkonen |45| who actually find C A 00 



0(e, W ,m 2 /i/ 2 )(Ci 
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Their estimate of the right- 



hand side of Eq. (|58ll does not contain our factor AA, 
but much more importantly they estimate /nl as if the 
right hand side were fully correlated with to conclude 
that the model can give /nl ~ 1- 

The reason of the discrepancy lies in the Eq. (71) in 
Ref. HH, which we write in the same form as our Eq. 



C A (t)-C 2 A (*i) 



eHM 2 Jti 



§HA~ 2 di{8&id i 8ax) 



12 In the proper treatment of the integral and its initial condition 
the factors e(£i) 2 , f/y(Cl) 2 , an d m &/H 2 (Ci) 2 cancel out since 



the evaluated quantity is in this case C^W 
conserved. This was not taken into account in Ref. 



C 2 A (*i 



and fi is 



Now we turn to the possibility that significant non- 
gaussianity could be generated during preheating. Pre- 
heating is the term used to describe the energy loss by 
scalar fields which might occur between the end of infla- 
tion and reheating |47t Esj . the latter being taken to cor- 
respond to the decay of individual particles which leads 
to more or less complete thermalisation of the Universe. 
Preheating typically produces marginally-relativistic par- 
ticles, which decay before reheating. 

It was suggested a long time ago 0, that preheat- 
ing might cause the cosmological curvature perturbation 
to vary at the level of first-order perturbation theory 
perhaps providing its main origin. More recently it has 
been suggested 0, |52 that preheating might cause 
the curvature perturbation to vary at second order, pro- 
viding the main source of its non-gaussianity. 

If the separate universe hypothesis is correct, a varia- 
tion of the curvature perturbation during preheating can 
occur only in models of preheating which contain a non- 
inflaton field that is light during inflation. This is not 
the case for the usual preheating models that were con- 
sidered in 0> EH E3 ' an d accordingly one does not ex- 
pect that significant non-gaussianity will be generated in 
those models 13 . This is not in conflict with the findings 



13 The preheating model considered in l4(i contains a field which 
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of pil I5H l52| because the curvature perturbation is not 
actually considered there. Instead the perturbation i/; MW 
in the longitudinal gauge is considered, which is only in- 
directly related to ( by Eqs. ©, O and 14 • We 
conjecture that non-gaussianity for the curvature pertur- 
bation on cosmological scales is not generated in the usual 
preheating models, but that instead the curvature per- 
turbation remains constant on cosmological scales. This 
should of course be checked, in the same spirit that the 
constancy of the curvature perturbation was checked at 
the first-order level p7j . 

The situation is different for preheating models which 
contain a non-inflaton field that is light during inflation. 
At least three types of models have been proposed with 
that feature Hi El El El- Except for [H only the 
magnitude of the curvature perturbation has been con- 
sidered, but in all three cases it might be that significant 
non-gaussianity is also generated. 



D. The curvaton scenario 

In the simplest version of the curvaton scenario 
the curvaton field a is ultra-light during inflation and has 
no significant evolution until it starts to oscillate during 
some radiation-dominated era. Until this oscillation gets 
under way, the curvature perturbation is supposed to be 
negligible (compared with its final observed value). The 
potential during the oscillation is taken to be quadratic, 
which will be a good approximation after a few Hubble 
times even if it fails initially. The curvature perturba- 
tion is generated during the oscillation, and is supposed 
to be conserved after the curvaton decays. Here we give 
a generally-valid formula for the non-gaussianity in the 
curvaton scenario, extending somewhat the earlier calcu- 
lations. 

The local energy density p a of the curvaton field is 
given by 



(63) 



where cr a (?7,x) represents the amplitude of the oscilla- 
tions and vn a is the effective mass. It is proportional to 
0,(77, x) -3 where a is the locally-defined scale factor. This 
means that the perturbation 8p a jp a is conserved if the 
slicing is chosen so that the expansion going from one 
slicing to the next is uniform [l5j. The flat sli cing cor- 
responding to ip MW = has this property 0, EI and 
accordingly 5p a is defined on that slicing. 

Assuming that the fractional perturbation is small 
(which we shall see is demanded by observation) it is 



may be heavy or light; we refer here to the part of the calculation 
that considers the former case. 
14 The slices of the longitudinal gauge are orthogonal to the threads 
of zero shear, and -0 MW on them is very different from the cur- 
vature perturbation 



given by 



Spa _ 2 Sa a 

Per <?a 



SOa 



(64) 



We first assume that cr(x) has no evolution between in- 
flation and the onset of oscillation. Then Sa a /a a will be 
equal to its value just after horizon exit, which we saw 
earlier will be practically gaussian. 

The total density perturbation is given by 



6_p 
P 



Spa 



(65) 



where Q a = p a j p oc a is the fraction of energy density 
contributed by the curvaton. Adopting the sudden-decay 
approximation, the constant curvature perturbation ob- 
taining after the curvaton decays is given by Eqs. 1|12|) 
and 114J1 . evaluated just before curvaton decay and with 
-0 = 0. In performing that calculation, the exact expres- 
sion Eq. 164f) can, without loss of generality, be identi- 
fied with the first-order part 5p ai /p aQ1 the second- and 
higher-order parts being set at zero. 

Adopting the first-order curvature perturbation in 
Eq. one finds E3 chi-squared non-gaussianity com- 
ing from the second term of Eq. (|64|l . 



NL 



5 

"47 ' 



with 



ip 



(66) 



(67) 



Ap r + 3p a 

evaluated just before decay (p r is the radiation den- 
sity) . Going to the second-order expression one finds jH^ 
additional chi-squared non-gaussianity. The final non- 
linearity parameter /nl = fm^ + 5/3 is given by 



5 5 5 



(68) 



If fl a -C 1 then /nl is strongly negative and the present 
bound on it requires CT > 0.01 (combined with the typ- 
ical value £ ~ 10~ 5 , this requires Sp a /p a « 1 as ad- 
vertised). If instead f2 ff = 1 to good accuracy, then 
/nl = +5/4. Either of these possibilities may be re- 
garded as generic whereas the intermediate possibility 
(|/nl| ~ f but /nl 7^ 5/4) requires a special value of fi CT 
just a bit less than 1. 

Finally, we consider the case that a evolves between 
horizon exit and the era when the sinusoidal oscillation 
begins. If a a (the amplitude of oscillation at the latter 
era) is some function g(a t ) of the value a few Hubble 
times after horizon exit, then 



5<r a = g'5a* + ^g"{5a*) 2 . 



(69) 



where the prime means derivative with respect to cr*. 
Repeating the above calculation one finds 



f 5 ^ 5 5 fu.99" 



(70) 
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The final term is the first-order result (given originally 
in the middle term is the second-order correction 

found in and the first term converts from f^if^ to 
/nl- 

E. The inhomogeneous reheating scenario 

The final scenario that has been suggested for the 
origin of the curvature perturbation is its generation 
during som e sp atially inhomogeneous reheating process 
EH EJL EM E3 ■ Before a reheating process the cosmic 
fluid is dominated by matter (non-relativistic particles, 
or small scalar field oscillations which are equivalent to 
particles) which then decay into thermalised radiation. 
At least one reheating process, presumably, has to occur 
to give the initial condition for Big Bang Nucleosynthesis, 
but there might be more than one. 

The inhomogeneous reheating scenario in its simplest 
form supposes that the curvature perturbation is negli- 
gible before the relevant reheating process, and constant 
afterwards. The inhomogeneous reheating corresponds 
to a spatially varying value (a perturbation) of the local 
Hubble parameter i? rc h(x) at the decay epoch (or equiv- 
alently of the local energy density), and this generates 
the final curvature perturbation. The perturbation in 
-ff rc h occurs presumably because it depends on the value 
of some non-inflaton 'modulon' field \ which was light or 
ultra-light during inflation. 

In contrast with the curvaton scenario, where the form 
p a can reasonably be taken as p a oc a 2 , the inhomoge- 
neous reheating scenario does not suggest any particular 
form for 7J re h(x)- Depending on the form, the inhomoge- 
neous reheating scenario presumably can produce a wide 
range of values for /nl- 

VIII. CONCLUSIONS 



with the current bound |/nl| < 10 2 seem relatively un- 
likely. In contrast, the simplest curvaton scenario can 
produce a strongly negative value (even violating the 
current bound). However, in the important special case 
where the curvaton dominates the energy density before 
it decays, it gives precisely /nl = +5/4. Finally, for the 
single-component inflation case, Maldacena's calculation 
combined with current constraints on the spectral tilt 
show that it has magnitude less than 1CU 2 . These result 
are summarised in the Table Q] 



TABLE I: Non-gaussianity according to different scenarios for 
the creation of the curvature perturbation. For the simplest 
curvaton scenario, /nl = +5/4 is a favoured value. 



Scenario 


I/nl| «1 


/nl ^ 


1 /nl < - 


1 /nl > 1 


Single-component 










inflation 


yes 


no 


no 


no 


Multi-component 










inflation 


likely 


possible 


possible 


possible 


Simplest 










curvaton scenario 


unlikely 


likely 


likely 


no 



In the near future, results from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) |6y| or elsewhere 
may detect a value |/nl| S> 1. If that does not hap- 
pen, then PLANCK 8] or a successor will either detect a 
value | /nl | ~ 1 , or place a bound | /nl 1^1- The precise 
level at which this will be possible has yet to be deter- 
mined because it would require a second-order calcula- 
tion of all relevant observational signatures. The exam- 
ple of the simplest curvaton scenario, where /nl = +5/4 
is a favoured value, shows that such a calculation and the 
eventual observations will be well worthwhile. 



We have examined a number of scenarios for the pro- 
duction of a non-gaussian primordial curvature perturba- 
tion, presenting the results with a unified notation. These 
are the single-component inflation, multi-component in- 
flation, preheating, curvaton, and inhomogeneous reheat- 
ing scenarios. Although the trispectrum may give a com- 
petitive observational signal 0j EH; we have focused 
only on the bispectrum which is characterised by the pa- 
rameter /nl- In all cases our treatment is based on exist- 
ing ones, though we do not always agree with the original 
authors. 

The preheating and inhomogeneous reheating scenar- 
ios cover a range of possibilities, which have not been 
fully explored but which can presumably allow a wide 
range for /nl- The same is true of multi-component in- 
flation, except that extremely large values comparable 
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